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Abstract. We study stability of isothermal two-component radiatively driven stellar winds against one-dimensional 
perturbations larger than the Sobolev length, and radiative-acoustic waves in such stellar winds. We perform 
linear perturbation analysis in comoving fluid-frames of individual components and obtain dispersion relation in 
the common fluid-frame. For high density winds the velocity difference between velocities of both components is 
relatively small and the wind is stable for radiative-acoustic waves discovered originally by Abbott, in accordance 
with the previous studies of the one-component wind. However, for such high density winds we found new types of 
waves including a special case of "frozen-in" wavy patterns. On the other hand, if the velocity difference between 
wind components is sufficiently large (for low density winds) then the multicomponent stellar wind is unstable 
even for large scale perturbations and ion runaway occurs. Thus, isothermal two-component stationary solutions 
of the radiatively line driven stellar wind with an abrupt lowering of the velocity gradient are unstable. 
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1. Introduction 

Since the foundation of the basic theory of the radiatively 
driven stellar wind the wind stability was one of the most 
fundamental questions to be solved. At the very beginning 
Lucy & Solomon ( 1970|) concluded that radiatively driven 
stellar winds are essentially unstable. Contrary to this as- 
sertion, Castor, Abbott &; Klein ( 1975|) described radia- 
tively driven stellar wind as smooth and stable steady- 
state outflow. This contradiction survived when Abbott 
( |1980| ) showed that the stellar wind described by CAK is 
stable. On the other hand, MacGregor et al. (1979) and 
Carlberg ( 1980 ) concluded that radiatively driven stellar 
winds are unstable. This paradox was solved by Owocki 
& Rybicki ( [1984 ) . These authors found a general relation 
which is valid for perturbation both smaller and larger 
than the Sobolev length, the so called "bridging relation" . 
The main result of the latter paper is that the flow of the 
line driven wind is stable for perturbations larger than the 
Sobolev length (the so-called large scale perturbations), 
yielding sta ble ra diative-acoustic waves, which were found 
by Abbott ( 1980D and unstable for perturbations smaller 
t han the Sobolev len gth, a s found by MacGregor et al. 
(|1979|) and Carlberg ( jl980| ). The theory of instabilities of 
radiatively driven stellar wind was further developed by 



Lucy fll984| ), Owocki & Rybicki (|l985| , |1986| , |l99l[ ), and 
extended to three-dimensional perturbations by Rybicki 
et al. ( |1990D . For an introduction to the problem of a sta- 
bility of a line driven wind see Rybicki ( |l987j ) and Owocki 
( 1992| ). The existence of instabilities in the wind is im- 
portant for an X-ray phenomenon because it is usually 
assumed that X-rays are generated by wind clumping or 
shocks (e.g Lucy & White |1980| , Lucy |1982| , Owocki & 
Cohen |l999| ). 

On the other hand, it is known that radiatively 
driven stellar winds have multicomponent nature (e.g., 



Springmann & Pauldrach 1992, Babel 1995, Porter 



Drew 1995). The stellar radiation is predominantly ab- 
sorbed by species (typically C, N, O, Fe, etc.) which have 
much lower density than the rest of the stellar wind, which 
is composed mainly of hydrogen and helium. However, 
multicomponent effects arc important only for low-density 
stellar winds. Recently, Krticka & Kubat fl2000| , ftOOla , 
|2001b[ hereafter KKO, KKI, KKII, respectively) computed 
models of isothermal two-component and non-isothermal 
three-component radiatively driven stellar winds. 
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Springmann & Pauldrach (1992) proposed that for low 
density stellar winds the absorbing component is not able 
to accelerate the non-absorbing component and that both 
components decouple. On the other hand, using a model of 
an isothermal two-component stellar wind KKO obtained a 
surprising result that the components do not decouple and 
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an unexpected decrease of the velocity gradient was found. 
This effect can be explained by the dependence of the 
radiative force on the velocity gradient. A question which 
naturally arises is the stability of such multicomponent 
flow. 

It is natural to expect that the conclusions about the 
stability of the one-component wind will be in principal 
also valid for the two component flows. Recently, Owocki 



& Puis (2002, hereafter OP) extended the general one- 
component stability analysis of Owocki & Rybicki ( 1984 ) 
for the case of a two-component isothermal wind and 
found that the two-component solution is unstable when 
t he fl ow is not well coupled. Here we extend Abbott's 
( 198C ) calculations to the case of the multicomponent flow 
and study how the overall picture of stable Abbott waves 
changes in the two-component isothermal stellar wind. 
The analysis presente d in t his paper is based on a part 
of a thesis of Krticka (|200l[). 



2. Time-dependent hydrodynamic equations for 
isothermal wind 

We assume an isothermal spherically symmetric wind 
consisting of two components, namely of passive (non- 
absorbing) hydrogen ions with mass equal to proton mass 
wip and charge equal to proton charge q p and of absorb- 
ing ions with mass Aim p and charge qi. Time-dependent 
radiatively driven stellar wind is then described by the set 
of hydrodynamic equations, namely with the continuity 
equations 



dpp 
dt 



l_d_ 

r 2 dr 



[r p p v ri 



dt r 2 dr ^ 1 ri ' 



= 0, 

= o, 



and with the equations of motion 



dv rp 

dv r i 

~dt 



dv rp 

v or 
dv ri 
or 



= -.9 



^rad 



Pp Pp or 

1 R ldpi 

- 9 -Rpi — 
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(la) 
(lb) 

(2a) 
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are velocities and 



In these equations v rp , p p , v ri , p\ 
densities of passive plasma and accelerated ions, respec- 
tively, p p , pi are partial gas pressures of each component 
(p p = a 2 p p , pi — afpi), isothermal sound velocities are 
a 2 v = kT/rrip and a 2 = kT/m u g = GM{1 - T)/r 2 is the 
gravitational acceleration (corrected for absorption by free 
electrons) acting on each component (G and 971 are grav- 
itational constant and stellar mass, respectively, T is the 
Eddington factor accounting for the absorption on free 
electrons) and g\ aA is the radiative acceleration acting on 
absorbing ions. We take the radiative acceleration in the 
form 



rad 



1 (J e L 
2); 47rr 2 c 



n e /W 



k 



2)i dv ri 
aeVthPi dr 



(3) 



with force multipliers k, a, S after Abbott ( 1982 ). Here / 
is the finite disk correction factor (Pauldrach et al. 1986, 



Friend & Abbott |1986| ), n e is the electron density (we 
set n e = rip), and W is stellar dilution factor. Here we 
have introduced the factor 2); (which is the ratio of the 
absorbing ions density to the passive plasma density in a 
stellar atmosphere) to account for radiative acceleration 
acting directly on ions (see KKO). 

Frictional force (per unit volume) R pl acting be- 
tween both components has following form (Springmann 
& Pauldrach |j~99^ ): 

-Rpi = -n p nikpiG(x pi ), (4) 

where n p and rt; are number densities of passive plasma 
and absorbing ions. The friction coefficient fc p i is given by 



47r In Aq 2 q 2 



k-n 



kT 



(5) 



where In A is the Coulomb logarithm, G{x) is the so-called 
Chandrasekhar function (see Springmann & Pauldrach 
1992| , KKO and Fig. for the run of this function). The 
argument x pl of the Chandrasekhar function in Eq. (||) is 
proportional to the ratio of the drift velocity 
to the thermal velocity v t h, namely 



1 r p 



^pi 



(6) 



where A pl — ApA-J (A p + Ai) is a reduced atomic mass. 



3. Radiative-acoustic waves 



Similarly to Abbott ( |1980D we keep the equations locally 
linear and we study waves in comoving fluid frames of 
individual components. Since the velocities of the compo- 
nents are different, we start with using two different fluid 
frames (one for each component). Thus, we use comoving 
fluid frames of non-absorbing and absorbing components 



r p =r' — Vp(r')t, 
ri = r' — Vi(r')t, 



(7a) 
(7b) 



instead of the frame of the static observer r' (note that 
both comoving fluid-frames are chosen to be local iner- 
tial fluid frames). Here both r p and n are radial coor- 
dinates. We assume that the wind is perturbed from its 
original steady-state and that the perturbed quantities do 
not change the density scale height. The perturbed quan- 
tities are denoted by 5p p , 6v p , 5p- u and Svi. 

The time-dependent continuity equations to the first 
order are 



dSpp 
dtp 
dSpi 
dU 



Po, P - 



Po,i- 



dSvp 
dr p 
dSvi 



dr. 



= 0, 



= 0, 



(8a) 
(8b) 



and two-component time-dependent linearized momentum 
equations are 

«p d5 P p G'(Av Q ) 



dSvr, 



dU 



Po, P dr p po.p G{Av ) 



(Svi ~ Sv p ) , (9a) 
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dSvi af dSpi 
dU p ,i dn 
R pi G'(Av a ] 



poa G(Av ] 



d v ,g iad 
(Sv : 



dSvi 
dri 



(9b) 



where the subscript denotes unperturbed quantity in 
the observer's frame, the velocity difference is Avq = 
vo,i - «o,p, G(Av ) = G(x ip ), G'(Av ) = dG(Av Q )/dAv 
and d v *g rad = dg\ ad /d (dv ,i/dri). We neglected gravity 
and density stratification in the momentum equation and 
assumed that the radiative force depends only on velocity 
gradient. In these equations we simply suppose Galilean 
transformation of coordinates for which t p = ti = t. We 
distinguish between t p and i; to emphasize difference in 
partial derivatives with respect to time. During calcula- 
tion of 8/ dtp one should keep r p constant and, similarly, 
for d/dti one should keep n constant. The Galilean trans- 
formation also implies that the velocity difference Avo and 
perturbations of velocities and densities Sp p , Spi, 5v p a Svi 
are the same in both inertial frames. We neglected density 
perturbations of the radiative and frictional forces. This 
is consistent with the calculations of Abbott ( 1980 ) where 
density perturbations of the radiative force have been also 
neglected. Justification of this neglect can be found, e.g., 
in the Appendix of OP. 

Taking a partial derivative of the equation (|8a| ) with 
respect to r p and substituting it to the partial derivative 
of the equation (|9a| ) with respect to t p gives 



d 2 Sv p 



2 d 2 Sv p 



R pi G'(Av ) 
Po, P G(Av ) 



dSvi dSv p 



(10a) 



Similarly, partial derivative of the equation (pq ) with re- 
spect to r; and subsequent substitution to the derivated 
equation (|^) with respect to ti gives 



d 2 Svj 



j d 2 Svi 
drf 



d v >g T 



;ad 



d 2 Svi 



R pi G'(Av ) 
po,i G(Avq) 



dndU 
dSvi dSv p 



(10b) 



We obtained two differential equations for velocity 
perturbations of both components. Unfortunately, these 
equations are written in different fluid frames. To pro- 
ceed further, we have to rewrite these equations in one 
common frame. We selected the fluid frame of accelerated 
ions. From the relations 



dr 2 drf ' 

d 2 

dt% ~ dtf ^dudn 



. 2 d 2 



(12c) 
(12d) 



Thus, we rewrite the Eq. (10a) as 



d 2 Sv D nA d 2 Sv v . nd 2 Sv r] 
2Aw 0T r-Tr Z + At'n 



dt 2 



dtidn 



drf 



2 d 2 5v p 
lp ~drf 



R pi G'(Avq) ( d d 
po, P G(Avo) \oti dn ' 



We assume a solution in the form of propagating waves, 
which in the reference frame of ions are 



Sv p = V p exp [i (uj p ti - K P n)] 
Svi — Vi exp [i (u>iti — Kin)] ■ 



(14a) 
(14b) 



The amplitudes V p , Vi are generally complex to account 
for phase shifts, and, since we are doing linear analysis, 
they do depend neither on n nor on ti. Similarly, u> p , uj{, 
k p , and Ki are independent of r; and t\. Substituting from 
( |l4| ) into the wave equations ([ll| |l0b| ) we obtain a system 
of equations 



(u p + Av q k p ) 2 V p - a p KpVp- 
.Rpi G'(Av Q ) 



Po, P G(Av ) 



(ujp + Av Kp) Vp 



3 i(w p t — KpTi) 



.Rpi G'(Av Q ) 
' p , P G(Av ) 



{uJi + AvoKOVie^*-^, (15a) 



[ufVi - afn?Vi + Kid v ,g rad UiVi- 



.R P i G'(Av Q ) 
Poa G{Av ) 

.Rpi G'(Av ) 

— — t 

p ,i G(Avq) 



"p'p 



We are looking for a non- trivial solution of Eqs. ( |15|)i i.e . 
a solution with Vi ^ 0, V p ^ for arbitrary t, n. Eq. (15a) 
can be rewritten as 

A pp V p e l ^ pt -^ n) = A pi Vie^ Uit ' K ' ri) 



or 



Vr, 



A p i 



exp [i (uii — uj p ) t — i (Ki — k p ) n] Vi. 



it follows that 

d _ d 
dr p dr-i ' 



n =rp - Av t p , 



— - — - Av — 

dt p dU dri ' 



(11a) 
(lib) 



(12a) 
(12b) 



Because V p , Vi, A pp and A p i do not depend on t and r;, 
the last equation can be fulfilled only if 



UJpt ~ Kp7*i = iO\t — KiH. 



(16) 



holds for any t, r; (the same conclusion can be obtained 
from Eq. ( |15b| )). Thus, wavenumbers and frequencies of 
both components are the same, 



Ki — K 



(17a) 
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(17b) 



Dividing Eqs. 



k p and lo — lo 



5|) by the exponential factor and denoting 
p we can rewrite the remaining system 



of equations as 



(oj + Avqk) 2 Ft, = al K 2 V n - 



■ p 



.i?ni G'(Av ) 



Po, P G(Av a ) 



(lo + Av q k) (V5 - V p ) , (18a) 



and we obtain the dispersion equation in the form of two- 
component Abbott waves 



(2 2 2 



)p ,p + (lo 2 + umd v >g 



rati 



2 2^ 
a ; K 



Poa 



0. 

(23) 



This equation corresponds to the one-component disper- 
sion relation and thus justifies the condition used by KKO 
to fix the mass-loss rate. Then solving the Eq. ( p3| ) for lo, 
the dispersion relation takes the form of 



uj 2 Vi = afn 2 Vi - Kd v ,g I3A uVi + 

.R P i G'(Av ) 
p ,i G(Avq) 



u(Vi-V p ). (18b) 



This system of equations can be rewritten in the ma- 
trix form as 



AV = 0, 



(19) 



where the vector V = (T/p,VS) T and individual elements 
of the matrix A are 



A p i 
An 



2^2 



2Av q hlo + AvqK 



2 2 

a p n 



iP p (lo + Av k) , 
iP p (ui + Av k) , 
iPiUj, 



2 2 2 
LO — CI: K 



where 



= R pi G'(Av ) 
p Po , P G(Av ) ' 



Pi = 



R pi G'(Avo) 



Po.i G(Avq) ' 
Equation ( p^| ) has a non-zero solution only if 
l|A||=0, 



(20a) 
(20b) 
(20c) 
(20d) 



(21a) 
(21b) 

(22) 



which is the dispersion relation. Generally, it has a com- 
plicated form and we will solve it numerically. However, in 
order to better understand the general dispersion relation 
we shall first find an analytical solution for some simpler 
specific cases. 

3.1. Abbott waves 

Let us assume that the velocity amplitudes of the com- 
ponents p and i are nearly equal (V\ w V p ) and that 
the phase velocity of the wave is much larger than the 
drift velocity (ui/k 3> Avq), i.e. the flow is well coupled. 



Then the terms containing Av in the Eq. (18a) can be ne- 



glected. These conditions mimic the one-component case, 
for which Abbott ( 198C| ) obtained stable radiative-acoustic 
waves, the so called Abbott waves. We may expect that for 
well coupled high density winds the two-component waves 
are similar to the Abbott waves. Indeed, summing disper- 
sion relations (fL8|) the imaginary frictional term vanishes, 



1 Poa 

2 PO.p 



'1 Am 

,2po, P 



dv>g r 



(24) 



where we neglected density of absorbing ions compared 
to the passive plasma density. Note that the Eq. ( p4| ) is 
the same as the Eq. (47) of Abbott ( |l980| ) , since we used 
a modified definition of the driving force Eq. (||) to ac- 
count for the fact that only the ionic gas is line driven. 
Apparently, lo is real and thus this mode is neither unsta- 
ble nor damped. There are two branches of lo correspond- 
ing to forward and backward waves, both depending on a 
wavenumber linearly. Note, however, that the general dis- 



persion relation derived by Owocki & Rybicki (1984) for 



the one-component flow allows for instabilities for short- 
wavelength perturbations. 

The point where the velocity of backward waves is 
equal to the wind velocity is called the critical point. The 
flow above the critical point cannot communicate with 
the wind base. Thus, early type stars are surrounded by 
the critical surface which separates two different domains 
of the stellar wind. Feldmeier & Shlosman ( 2000 ) aptly 
compare this situation to the cosmic censorship hypoth- 
esis. The result that the multicomponent nature of the 
wind does not alter the Abbott waves is important be- 
cause these waves determine the mass-loss of the CAK 



wind (Feldmeier & Shlosman 2001) 



Finally, in the case when the gas-pressure term a p can 
be neglected, the dispersion relation (|24|) corresponds to 
the lo+ mode of OP, namely 



LO + 



— k o v ig 



(25) 



We denoted this mode in accordance with OP as 
lo + although this mode is upstream. Note that 
Eq.(24) has another solution for a p = 0, namely 

lo 



0. 



3.2. Purely ionic Abbott waves 

In t he ca se when |V1| ^ \V P \ we can obtain from the 
Eq. ( |l8b| ) the dispersion relation in the form of Abbott 
waves in absorbing ions 



2 2 2 

LO = CI : K 



lok d v *g 



rad 



' Po,i G(Av ) 



(26) 
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Fig. 

Eq. 



0.25 



0.2 



0.15 



0.1 



0.05 



G'>0 
stable 



G'<0 
unstable 



1. The run of Chandrasekhar function (note, that 
yields x\ v ~ Avq). When the flow is well coupled, 



Xi p < 0.97, G'(Avo) > and the wind is stable. When 
the drift velocity is large, x lv > 0.97, G'(Avo) < and 
the wind is unstable. Note that the point a^p w 0.97 cor- 
responds to the maximum of G. 



not the case, because the imaginary term in the Eq. ( |l8b| ) 
is larger than in Eq. (18a) and thus Eq. (18b) does not 



allow |Vp| ^> |Vi|. Therefore there are no acoustic coun- 
terparts of such waves in a non-absorbing component, i.e. 
there are no passive plasma waves for which the condition 
|Vp| > \Vi\ is valid. 

3.3. Zero factional force 

To complete the list of simplified cases, we have also to 
mention the hypothetical case when interaction between 
the components vanishes. If both flow components do not 
influence each other, i.e. if the frictional force is zero, 
then the system of equations ([l^) does not implicate that 
the frequencies and wavenumbers of both components are 
equal. Instead of the system of dispersion relations (18) 
we obtain two independent relations for each wind com- 
ponent 



(w p + Av k p ) - a p K p = 0, 
md^c^Ui = 0. 



2 2 2 

UJ-. — a-. K- 



(28a) 
(28b) 



We will study such waves in the case when the pres- 
sure term a?K 2 in the dispersion equation ( p6| ) can be ne- 
glected. In this case, the dispersion relation has a form, 
which corresponds to the w_ mode of OP, 



The dispersion relation of the nonabsorbing component 



ad 



; R pi G'(Avq) 
Vo.i G(Ai> ) ' 



(27) 



The real part of cj_ depends on n linearly, whereas the 
imaginary part does not depend on k. Due to the pres- 
ence of the imaginary term in Eq. (|27]) such waves are 
damped in the case when G'(Avq) > 0. For clearness we 
plotted the run of the Chandrasekhar function in Fig. [j] 
Thus, if the wind is coupled [when the drift velocity is 
lower than the thermal velocity or, more precisely, when 
the Chandrasekhar function is a rising function of Avq 
(G'(Avq) > 0)], the wind is stable for this type of waves. 
On the other hand, when the Chandrasekhar function 
G(x) is decreasing (G'(Avq) < 0), the wind is not sta- 
ble for the above mentioned perturbations, and this mode 
leads to an ion runaway. 

Physical reasons for such instability are straightfor- 
ward. If the Chandrasekhar function is before its maxi- 
mum (as function of Avq), then the increase of the veloc- 
ity difference between velocities of ions and passive plasma 
enhances the frictional force, which finally tends to lower 
the velocity difference yielding a stable flow. In the oppo- 
site case (if the velocity difference is larger than that cor- 
responding to the maximum of Chandrasekhar function) 
the increase of the velocity difference lowers the frictional 
force which allows for additional increase of the velocity 
difference. Such two-component flow is clearly unstable. 
This is the ef fect o f the ion runaway (see also Springmann 
& Pauldrach [l992] OP). 



At a first glance there might exist similar acoustic 
waves in the non-absorbing component. However, this is 



(28a) describes ordinary isothermal sound waves and the 
dispersion relation of the ionic component (28b) describes 
stable Abbott waves of absorbing ions. 

Note that a similar result of independent waves can be 
obtained for the maximum of Chandrasekhar function, for 



which G'(Avq) — and the interaction terms in Eqs. (18) 
vanish as for the case of zero frictional force. 



4. Numerical results 



The simplified calculations presented in Sections 3.1 and 



^2] can help us to better understand the behavior of in- 
dividual branches of the general dispersion relation (^2|) . 
These calculations have been done numerically. To study 
the individual branches of the dispersion relation we 
solved numerically th e Eq. ( ^2|) using the procedure CP0LY 
of Jenkins & Traub ( |l970|) . 



Moreover, solving Eq. (18a) or Eq. ( 18b ) for given to 
and k we obtain the relation between wind amplitudes 



Vs = - 



A ■ p 



or, equivalently, 



Vi = --^ 



(29a) 



(29b) 



4.1. Stable wind for ionic Abbott waves 

First, we shall study the case of the star e Ori described 
in KK0. This star has a relatively dense wind, so the drift 
velocity between both components is low compared to the 
thermal speed of hydrogen and thus this wind should be 
stable for ionic A bbo tt waves because G'(Avq) > in this 
case (cf. Section 3^). To be specific, we shall study the 
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Fig. 2. Individual branches of the dispersion relation for 
the point where v r ~ 0.3 Voo in the dense wind for a spe- 
cific case of e Ori. Upper panel: Real roots of u>. Solid line 
represents a forward wave, dashed ones backward waves. 
Note, that the real part of branch denoted as (4) changes 
its signature and for k = kq is the real part of frequency 
u> zero. Lower panel: Imaginary roots of w, all waves are 
damped. As in the upper panel, solid and dashed lines de- 
note forward and backward waves, respectively. Real and 
imaginary parts of individual roots are denoted using the 
same numbers. 



dispersion relation in the point where v r rj 0.3 Voq. This 
choice does not influence the overall picture of our results, 
however. 

Resulting dispersion relations are displayed in the 
Fig. ^. Real parts of branches show linear dispersion re- 
lations as was predicted in Eqs. (|24|) and (p7|). The real 
branches of Abbott waves (two middle straight lines in 
both panels of Fig. |^ denoted by (2) and (3)) in the two- 
component case are essentially the same as in the one- 
component case (cf. Abbott 1980] ). From the correspond- 
ing imaginary parts of these two-component Abbott waves 
follows that these waves are damped only marginally. If 



we extend the calculations to n > 10 cm we could 
in principle obtain larger damping, but these values are 
beyond the region of validity of the assumption of per- 



turbations larger than the Sobolev length. We also calcu- 
lated the relation between wave amplitudes Vj and V p o f 
absorbing and nonabsorbing components using Eq. ( |29a| ). 
This calculation confirmed the assumption used in Sect. 
|3.l| for calculation of two-component Abbott waves that 
wave amplitudes for this waves are nearly the same. 

Largest imaginary branch corresponding to the ionic 
Abbott waves (upper straight line in both panels of Fig. |2[) 
does not depend on k (see Eq. (]27|)). Clearly, accord- 
ing to previous results, such waves are heavily damped. 
Moreover, the calculation of the relation between V\ and 
V v confirmed that |VJ| 3> |Vp| for these waves (cf. Sect. 

m 

However, a new type of slow waves appeared. They 
are described by the lower curves in both panels of Fig. || 
(denoted as (4)). The real part of the dispersion relations 
shows almost linear dependence with the exception of the 
region where it passes through the value of 9?(u>) = 0. 
Consequently, this solution corresponds to both forward 
and backward waves, which are stable and only marginally 
damped. The case of 3?(w) = deserves special attention. 
It corresponds to a static wavy structure in the comoving 
frame, so in the observer frame this structure resembles 
almost stable outflowing " frozen-in" wavy patterns of the 
characteristic size of n^ 1 w 10 7 cm for this specific case 
(note that kq is wavenumber for which 5R(w) = 0). 

In addition, calculations showed that the value of 
wavenumber kq (for which 5R(cj) = 0) depends on the dis- 
tance from the star. At the base of the wind the value of kq 
is lower, kq ~ 10 _6 cm _1 implying the possible character- 
istic pattern size of the order 10 6 cm whereas in the outer 
parts of the wind kq rj 10~ 10 cm -1 yielding the character- 
istic pattern size of the order 10 10 cm. 

For this new type of waves we can obtain approxi- 
mate analytical dispersion relation. Calculations showed 
that due to the low value of ui for these waves all terms 
in the dispersion relation (|22] ) can be neglected except 
constant terms and terms linear in ui. Further neglect of 
all terms which do not significantly influence this type of 
waves leads to the dispersion relation in the form of 



-a 2 d v ,g™ d n + ialPi) + n 2 a 2 a 2 + iAv P p afK 



0. 

(30) 



The condition JJ(a;) = can be now written as 

K 2 a 2 p d v ,g lad - AvoPpPi a (31) 

from which we can obtain the equation for kq in the form 
of 



1 AvoPpPi 



(32) 



However, we must keep in mind that our analysis was 
only linear, taking into account the nonlinear effects may 
change this promising picture, similarly to the case of 
one-component Abbott waves which become unstable if 
the s econd order effects are taken into account (Feldmeier 
1998b . 
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Fig. 3. Dispersion relation for the point r = 2.4 i?* in the 
wind of a B5 star with q\ = 4.8g p . Upper panel: Real roots 
of u). Solid line represents a forward wave, dashed ones 
backward waves. The relation approximately correspond- 
ing to the dispersion relation u> = —Avqk (see Eq. (|33|)) 
is denoted. Lower panel: Imaginary roots of ui, all waves 
lead to an instability. As in the upper panel, solid and 
dashed lines denote forward and backward waves, respec- 
tively. Real and imaginary parts of individual roots are 
denoted using the same numbers. 

4.2. Unstable wind for ionic Abbott waves 

Second we study a model of a B5 star with artificially 
enhanced effect of friction (g; = 4.8g p , see KKO) where 
the multicomponent nature of the flow plays an impor- 
tant role and leads to lower outflow velocity than in a 
one-component case. Here we study the dispersion rela- 
tion at the point r — 2.4 i?* where the drift velocity be- 
tween both components exceeds the value corresponding 
to the maximum of Chandrasekhar function and, there- 
fore, G'(Av ) < (see Fig. [TJ) _ 

The dispersion relations displayed in the Fig. || sub- 
stantially differs from the previous case. Consistently with 
simplified considerations in the Sect. ^2, the imaginary 
parts of all roots are negative and therefore the flow is 
unstable. The branch with the largest absolute value of 
both real and imaginary parts corresponds to the ionic 



Abbott waves. According to the Eq. (p7|), the real part of 
this branch depends on k linearly whereas the imaginary 
part does not depend on k. 

A natural question arises: Where have the two- 
component Abbott waves disappeared? However, the sim- 
plified calculations yielding the two-component Abbott 
waves (see Sect. 3.1) are not valid in this case, because 
the velocity difference Ai>o cannot be simply neglected. 
On the other hand, the dispersion relation ( |22] ) can be ap- 
proximately fulfilled if the term uj + Avqk vanishes. This 
conclusion can be simply justified because the dispersion 
relation (E2) can be rewritten as 



A, 



A- — A ■ A- — n 

Avqk — then the ^4 pi 
tically and A pp vanishes if the pressure term a p K" is neg- 
ligible (see Eqs. (20a, 20b| )). The linear dispersion relation 



and when 



term vanishes iden- 

2 J2 



uj = — Avqk (33) 
is visible in the Fig. || (note that Avo « 10 6 cms -1 ). 

4.3. On the stability of the B5 star wind model 



Nevertheless, the analysis presented in the Section 4.2 
does not mean that the whole wind is unstable for the 
above mentioned two-component instability. Since stabil- 
ity depends on the sign of G'(Avq) which varies through- 
out the flow, there are regions where the wind is stable 
and regions where the wind is unstable. This situation is 
depicted in the Figure ^. 

Near the stellar surface the wind is well coupled 
(uj/k S> Avq), the Chandrasekhar function is below 
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Fig. 4. The stability of the wind of a B5 star with q\ = 
i.8q p . At the base of the wind x lp < 0.97 (note, that 
the value Xi p w 0.97 approximately corresponds to the 
maximum of Chandrasekhar function, see Fig. [l]) and the 
wind is stable. In the outer parts the density decreases, 
Chandrasekhar function reaches its maximum and Xi p > 
0.97, the wind is not stable. In the outermost parts of the 
wind is x lp < 0.97 again and the wind is stable. 
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its maximum (as a function of the velocity difference), 
G'(Avo) > and the wind is stable (see Eq. @). This 
conclusion was justified by numerical calculations at the 
point where v rp « 4a p (obtained dispersion relation re- 
sembles that in Fig. ^). 

Downstream is the wind accelerated, its density de- 
creases and the velocity difference Avq between both com- 
ponents increases. However, such low density wind reaches 
the point where the maximum of the Chandrasekhar func- 
tion is reached and the absorbing component is not able to 
accelerate the passive component sufficiently. Due to the 
functional dependence of the radiative force the velocity 
gradient of both components decreases (see KKO and the 
discussion in OP), and, in addition, G'( Avq) < 0. As has 
already been shown in the Section 12 (see Fig. ||), the 
flow is unstable in this region. 

In the outermost parts of the wind the relation 
G'(Avq) > holds again and clearly, the wind is stable 
there. However, because the wind is unstable upstream 
where G'(Avo) < 0, the instabilities from that unstable 
region may disseminate and influence the stability of the 
outermost parts of the wind. Such effects would be studied 
using hydrodynamical simulations. 

Summarizing, there are three regions in the wind of 
this B5 star, namely the innermost and outermost parts 
of the wind are stable against perturbations larger than 
the Sobolev length, whilst the "middle" part of the wind 
is unstable (see Fig. |^). 

The new slow waves that appeared in the stable solu- 



tion for a dense wind of an O star (cf. Section 4.1) are 
also present in the stable parts of the wind of a B5 star. 
For the inner region of stability these waves may be both 
forward and backward with a special case of "frozen-in" 
waves for k = no where 5R(w) = 0. Similarly to the case of 
an O star, the value of kq decreases for increasing radii. 
These waves are also present in the outer stability region, 
but in this case these waves are purely forward. 

4.4. The stability of a B5 star model with normal 
friction 

Finally, we studied the stability of a wind model of B5 
star with q\ = 2q p , which is a bit closer to real winds of 
B5 stars, at the point where v ri « 290 kms -1 . At this 
specific point is x- lp * 1.4 and thus, the flow is unstable 
there. The results of this calculation are similar to those 
displayed in the Fig. |[ All imaginary roots are negative, 
the flow is unstable at this selected point. The analysis 
presented in the preceding section is valid also for this 
case. 

5. Conclusions 

We showed that two-component isothermal radiatively 
driven stellar wind is unstable in the case when friction 
affects the overall structure of the wind, i.e. when the drift 
velocity between both components is sufficiently large. 
Strictly speaking, the wind is unstable if the argument of 



Chandrasekhar function is larger than the value for which 
the maximum of Chandrasekhar function is reached. For 
this case the so-called ion runaway instability occurs. In 
the opposite case, when the wind is well coupled, the 
Abbott waves (for large scale perturbations) in the wind 
are stable. Putting these two cases together, stationary 
wind solutions obtained by KKO for the case when the drift 
velocity increases to such values that the Chandrasekhar 
function passes through the point with its maximum value 
(i.e. it is decreasing and its derivative is negative) are not 
stable. Note that the region of instability falls within the 
region of the abrupt decrease of the velocity gradient in 
the solution found by KKO. 

In the real case the obtained large growth rate of 
the instability will be probably reduced by the en- 
ergy dissipation via frictional heating. Nevertheless, the 
inclusion of frictional heating will probably not alter 
the presence and the overall picture of instability for 
high drift speeds. However, it is not clear whether the 
ionic component escapes the star separately (as proposed 
Springmann & Pauldrach 1992 ), because the ionic Abbott 
waves have only modest spatial growth rate (see OP). 
Hydrodynamical simulations are necessary to resolve this 
problem. However, it is questionable whether the Sobolev 
approximation may be used for such calculations. The 
Sobolev approximation was at the edge of its validity dur- 
ing the calculations of KKO. If the ionic component leaves 
the star separately, then the Sobolev approximation can 
be used thanks to a large velocity gradient of absorbing 
ions after the decoupling (see KKI). On the other hand, 
if the decoupling process is more complex with modest 
velocity gradients, we may come beyond the region of va- 
lidity of the Sobolev approximation. 

More complete dispersion relations for two-component 
flow should be derived with the inclusion of the short- 
wavelength perturbations. In the one-component case 
these perturbations lead to the well-kn own r adiatively 
drive n wind instability ( MacG regor et al. 1979 , Carlberg 



1980| , Owocki & Rybicki |l984[ ) . When such perturbations 



on a scale below the Sobolev length are included, then the 
onset of a two-component instability occurs even below 
the maximum of Chandrasekhar function (see OP). 

The existence of the two-component instability could 
also lead to frictional heating of the wind up to the tem- 
peratures of the order of 10 6 K (cf. KKII). Such heating 
could explain enhanced X-ray activity of many B-stars. 
Similar mechanism of X-ray generation was also proposed 
by Porter & Drew ( 1995| ) and OP. Unfortunately, in the 
latest numerical models of a three-component nonisothcr- 
mal wind of B stars the velocity difference is too low to al- 
low for such ion runaway instability (except for extremely 
low density winds). On the other hand these results are 
based on a bit artificial dependence of the radiative force 
on the temperature and a more advanced calculations can 
alter this result. 

In a case of a wind where two component effects are 
not important (an O star wind), the two-component sta- 
bility analysis enabled us to find more types of waves than 
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the classical one-component analysis of Abbott ( |1980| ) . In 
addition to the original stable Abbott waves, there ex- 
ist heavily damped ionic waves and very slow waves with 
very weak damping. The latter waves may move in both 
directions with respect of the wind or they may be static 
in the comoving fluid frame. Such slow waves resemble 
"frozen- in" wavy pattern of the dimension of approxi- 
mately ~ 10 7 cm. 
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